Dynamical Mean-Field Theory for Quantum Chemistry 
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The dynamical mean-field concept of approximating an unsolvable many-body problem in terms 
of the solution of an auxiliary quantum impurity problem, introduced to study bulk materials with a 
continuous energy spectrum, is here extended to molecules, i.e. finite systems with a discrete energy 
spectrum. The application to small clusters of hydrogen atoms yields ground state energies which 
are competitive with leading quantum chemical approaches at intermediate and large interatomic 
distances as well as good approximations to the excitation spectrum. 
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Determining the ground state energy and dynamical 
properties of a system of N interacting electrons is a fun- 
damental unsolved problem. The properties of systems 
with "strong correlations," those dominated by local in- 
teraction effects, have been proven particularly difficult 
to treat. Over the years many approaches have been de- 
veloped fH4l5j|. but no one method has emerged as gen- 
erally applicable. In particular, for all but the smallest 
molecules the treatment of correlation effects in quantum 
chemistry (for example, bond breaking, or the energetics, 
dynamics, and magnetic properties of transition metal 
clusters) remains a frontier research area. 

In recent years the theoretical study of strongly cor- 
related condensed matter systems has been revolution- 
ized by the development of dynamical mean-field theory 
fDMFT), first in its single-site 0, E3] and then in its 
cluster [18l42fj | forms. DMFT is a Green's function based 
method, in contrast to many of the quantum chemical 
methods which are wave function based, and is often pre- 
sented in terms of an impurity self-consistently coupled to 
a noninteracting bath of states. To date, the main appli- 
cations 0,0, EH of DMFT have been to extended (typ- 
ically periodic although an interesting recent application 
to a nanoscopic conductor should be noted systems, 
characterized by a continuous density of states. Here, we 
show that DMFT can be used for finite systems, for which 
the "bath" is characterized by a discrete (even small) 
density of states. Application to a benchmark quantum 
chemical system (H„, the n- hydrogen molecule in vari- 
ous configurations) suggests the method may be useful 
for treating the strong correlation problems of quantum 
chemistry. 

In the same way that density functional theory is 
derived from the Hohenberg-Kohn density functional 
DMFT may be derived from the Luttinger-Ward 
functional $lw [HI, which is a functional of the elec- 
tron Green's function G, the interparticle interaction 
I and an external potential V(r). G is defined in 
terms of the operator i]) a (r,t) which annihilates an elec- 



tron at position r and time t as G aa t{v,v' ,t — t') = 
T t Vct (r, t) , i\)\, {r', t') ^. The appropriate G for a given 

external potential V satisfies S&ujv/SG = 0. 

In analogy to density functional theory, $lw may be 
written as the sum of two terms: a "universal" term <f>L\v v 
defined in terms of the sum of all vacuum to vacuum 
Feynman diagrams which has explicit dependence only 
on G and I (not on V) and a material-specific term which 
depends explicitly on V and G but not I. DMFT is an 
exact extrcmization of an approximation to $lw V j m the 
same way that practical implementations of density func- 
tional theory are exact extremizations of approximations 
(such as the local density approximation) to the exact 
density functional. 

In practice most correlated electron calculations pro- 
ceed by reducing the full problem (which may be viewed 
as a matrix in the space spanned by the complete set 
of states </> a (r)) to a "correlated" subspace spanned by 
a set of correlated states, (/>™ rr (r). One defines a corre- 
lated problem by retaining only the matrix elements of 
G and I within the correlated subspace and writing a 
Luttinger-Ward functional for the correlated degrees of 
freedom as 

*lw = - Trln [Gq 'Gcorr] + Tr [G^G^,.] (1) 



Go = (idtl — H£| rr ) _ is the noninteracting Green's 
function defined in the usual way in terms of the <p a and 
restricted to the subspace of correlated states while &™" v 
is formally defined as the sum of all vacuum to vacuum 
diagrams (with appropriate symmetry factors) involving 
G corr and interactions I corr . 

Within the correlated subspace we define the dynam- 



ical self-energy S = G^ 



G: 



If the correlated 



subspace contains M states (summed over atoms and 
orbitals), £ may be represented as M(M + l)/2 func- 
tions of frequency. The DMFT method approximates the 
self-energy as a sum of a much smaller number of func- 
tions. Different versions exist [3, 18, 19j], corresponding 
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to different approximations to S. Each approximation 
to S implies an approximation to the Luttingcr-Ward 
functional, which is such that the extremization may be 
carried out by solving a quantum impurity model with 
parameters fixed by a self-consistency condition. For 
specifics, see Ref. [13, H3|- The results presented here 
were obtained with the cellular DMFT (CDMFT) ver- 
sion , which has a real space interpretation naturally 
suited to the nontranslationally invariant problem posed 
by a molecule. 

In CDMFT one divides the full set of states into P 
cells, labeled by cell index J = 1, ...,P, with each cell 
containing some small number N of orbitals, labeled by 
an orbital index a = 1,...,N. (Note that each cell may 
contain more than one atom.) The total number of cor- 
related states is M = PN . The approximation is to re- 
tain only those components £( Jiai )w2a 2 ) of S for which 
Ji = Ji (i.e., the two orbitals are in the same cell) and 
correspondingly only the intracell terms in I corr . Defin- 
ing the resulting self-energy as Xcdmft we obtain an 
approximation to the Green's function: 
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An approximate Luttingcr-Ward functional is con- 



structed by using G a 



in the Tr and Trln terms in 



Eq. (JTJ and approximating the universal part as the sum 
of P iV-orbital quantum impurity models (different in 
each cell in nonperiodic cases) defined in terms of a local 
Green's function and the intracell portions of the inter- 
action I corr . The Sq DMFT are found from the solution of 
the quantum impurity model for cell J. The stationarity 
condition on the approximate Luttingcr-Ward functional 
is that the quantum impurity model Green's function for 
cell J equals the J — J component of G approx [Eq. fl2J] ; 
enforcing this condition fixes the parameters of the quan- 
tum impurity model [171 ]. The procedure becomes exact 
as cluster size N approaches system size M and provides 
a reasonable approximation for the small values of N 
which are computationally accessible [24j. Nothing in 
this derivation requires that Go arise from a system with 
P = oo and a continuous density of states. 

We apply the CDMFT approximation to a standard 
quantum chemical reference system, the H„ molecule 
consisting of n hydrogen atoms which was recently stud- 
ied by Tsuchimochi and Scuseria [Toj using a constrained- 
pairing mean-field theory (CPMFT). We present here re- 
sults for H„ chains and rings with n = 6 and n = 50, 
as well as the H4 tetrahedron. By chain we mean a lin- 
ear arrangement of atoms with open boundary conditions 
and interatomic spacing i?, and by ring we mean a circu- 
lar arrangement with chord distance R between nearest 
neighbor atoms. We follow Ref. and define the cor- 
related subspace as the set of hydrogen Is orbitals (with 
both spin directions) centered on each of the hydrogen 
atom positions R a : <\) a ^ a = ipi s (R — R a )- The interaction 



-0.48 



-0.5 



-0.52 



I" 



/ // 



H fi chain 



-0.54 



-0.56 - 



Exact(AII) 



UHF 

Exact(U) - *- 
ExacttU.vJ ■■•»■- 
DMFT(U) — «■- - 
CDMFTfU.vJ - .- 



2 3 4 

R H . H (bohr) 



-0.48 




-0.52 



-0.54 - 
-0.56 - ^ 



H 6 ring 



Exact(AII) — 
UHF -~- 
Exact(U) - ■*- 
ExactfU.V,) --T-- 
DMFT(U) — - 
CDMFT(U) - •- 
CDMFTfU.v, --. - 



R H - H (bohr) 

FIG. 1: Energy of Hg chain (upper panel) and ring (lower 
panel) as a function of interatomic spacing R calculated with 
different methods as indicated. 



porr j g obtained from the appropriately antisymmetrized 
Coulomb integrals, which we compute using the standard 
minimal STO-6G basis [HI . The long range of the inter- 
action means that all sites are coupled. In practice one 
can treat accurately only a subset of the interactions; the 
others must be treated more approximately. We present 
energies in standard atomic units (a.u.); in the STO-6G 
basis the energy of the isolated H atom is -0.471 a.u. We 
use unrestricted Hartrce-Fock (UHF); other choices (e.g., 
density functional theory) are also possible [2l| . 

To determine which interactions must be retained we 
observe that for n = 6 the problem posed by the cor- 
related (Is) subspace with the complete interaction can 
be diagonalized exactly. The upper panel of Fig. [TJ com- 
pares the resulting energies to those obtained by treat- 
ing some parts of the interaction via UHF and the other 
parts exactly. We focus on the on-site term, convention- 
ally denoted U = I^a'a a an d the first neighbor terms 

I™a,b,b> I C a°Aa,b with b a nearest neighbor of site a, 
which we lump together into a term V±. We define two 
models, "Exact(£/)," where we treat the U term exactly 
and all other interactions by UHF, and "Exact (t/, Vi)" 
where we treat the U and all V\ terms exactly and all 
other interactions by UHF. From Fig. Q]we see that the 
Exact (Z7, Vi) approximation is much closer to the exact 
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energy than the Exact(J7) approximation. We conclude 
that it is important to use a method which incorporates 
the on-site and first neighbor interactions, while the re- 
maining terms may be treated approximately. 

Each site has two (counting spin degeneracy) orbitals, 
so as defined above the CDMFT method with N = 2 
means the cluster contains one atom. At this level we 
have an approximation to the Exact([/) Hamiltonian 
which we refer to as "DMFT([/)." The CDMFT method 
with N = 4 corresponds to a two-atom cluster and consti- 
tutes an approximation to the Exact([/, Vi) Hamiltonian 
which we refer to as "CDMFT([7, Vi)." We have also used 
CDMFT as an approximation to the Exact (t/) Hamilto- 
nian. We refer to this approximation as "CDMFT(Z7)-" 

We have used these methods to calculate the ground 
state energies of the Hq chain and ring. In most of our cal- 
culations we have used the "exact diagonalization" (ED) 
method ptjj to solve the quantum impurity model. Up 
to 10 bath sites were used in the CDMFT (1/, V\) cal- 
culation of the H50 chain; convergence with a number 
of bath sites was verified. In a few cases we also veri- 
fied that a continuous-time quantum Monte Carlo (CT- 
QMC) method [23| gives identical results. The two panels 
of Fig.Q]sIiow that each approximation produces a result 
which lies somewhat above the exact result for the Hamil- 
tonian which it approximates, with CDMFT providing a 
better approximation than DMFT. Also, although it is 
difficult to perceive in the figures, the DMFT(J7) equa- 
tions have a phase transition to an antiferromagnetic 
state at R <~ 4 bohr while the Exact (U,Vx) does not 
have any phase transition. Both DMFT approximations 
locate the minimum in the E(R) curve at essentially the 
exact position (~ 1.75 bohr), unlike the UHF approxima- 
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FIG. 2: Main panel: Energy of H50 chain as a function 
of interatomic separation R calculated by the methods of 
this Letter and compared to digitization of data presented 
in Ref. [l(J for coupled cluster methods (CCSD), DMRG 
(Ref. [H) and other wave function based methods (CPMFT). 
Inset: Relative error of CPMFT and DMFT methods, defined 

as (-E mo thod — -EdMR.g)/(-BdMRG — -^isolated atom). 



tion (~ 1.7 bohr). A recent Hubbard model study [24j 
found that going from single-site to two-site clusters im- 
proved the energy substantially; larger cluster sizes con- 
verged slowly to the exact result and the differences be- 
tween two-site and exact results are less than the errors 
involved in constructing the approximate Hamiltonians, 
and will not be of interest here. 

The main panel and inset of Fig. [2] compare the ener- 
gies of a 50 atom hydrogen chain obtained by our single- 
site and two-site cluster DMFT methods to the ener- 
gies obtained from other approximate methods [To| and 
from the density matrix renormalization group (DMRG), 
which is believed to provide a numerically exact solution 
to the problem. The solid line is a coupled cluster cal- 
culation also reported in Ref. [l(|; this weak coupling 
method fails to converge as the separation increases be- 
yond R > 2.1 bohr. The dotted and dashed lines la- 
beled CPMFT are obtained from a variational wave func- 
tion method The inset, which shows errors relative 
to DMRG, demonstrates that even the relatively primi- 
tive DMFT methods used here are gencrically accurate 
over the whole intermediate to strongly correlated regime 
(R > 1.5 bohr). In the weakly correlated {R < 1.5 bohr) 
regime, other methods (for example coupled cluster) are 
preferred, but these methods fail in the strongly corre- 
lated (here R > 2.1 bohr) regime. 

Figure [3] shows energies obtained for the tetrahedrally 
coordinated H4 molecule, which is methodologically chal- 
lenging because it is only weakly bound. Here UHF fails 
qualitatively, predicting that the molecule is not bound 
at all. DMFT(Z7) is qualitatively worse, even though the 
intcrsite terms in the Hamiltonian make only a small con- 
tribution to the ground state energy, as is seen from the 
close correspondence of the Exact (U) and exact energy 
traces. Remarkably, the CDMFT(J7, Vi) trace [applied, 
of course, to the Exact (£/, V\) model] produces a reason- 
able approximation to the energy, lying much closer to 
the exact curve and, in particular, predicting a minimum 
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FIG. 3: Energy curves of tetrahedrally configured H4 calcu- 
lated by different methods as indicated. 
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FIG. 4: Main panel: Electron addition spectra for He ring 
calculated from exact diagonalization of the full Hamiltonian 
in the correlated subspace and from the CDMFT([/, Vi) equa- 
tions using both ED and CT-QMC solvers. Inset: Full fre- 
quency range. 



provide a new context in which to examine basic theo- 
retical issues including the "double counting correction" 
needed when combining dynamical mean-field and den- 
sity functional theory, as well as the possibility [2l| of us- 
ing dynamical mean-field self-consistency ideas to more 
systematically define the correlated subspace. Research 
in all of these directions is in progress. 
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(although not quite at the correct location) even though 
the UHF curve does not. This illustrates that appro- 
priately chosen cluster methods can capture even quite 
subtle behavior. 

Figure S] presents the electron spectral function (many- 
body density of states) projected onto one site of an Hg 
ring. The spectra are discrete, and the exact solution 
and the ED curves have been artificially broadened. For 
the CT-QMC curve we have used the maximum entropy 
methods of Ref. [28| without any additional broadening. 
We sec that CDMFT([7, Vi) (with cither solver) provides 
a good approximation to the exact curve, reproducing the 
gap in the excitation spectrum and the basic structure of 
the electron addition and removal spectra. 

In summary, we have shown that in the intermediate to 
strong correlation regime, appropriately chosen DMFT 
methods give results which are superior to other approx- 
imate methods. This, and the observation that the com- 
putational cost scales linearly with the system size (with 
a prefactor which depends strongly on the DMFT clus- 
ter), motivates a broader exploration of DMFT based 
methods in the quantum chemical context. A more sys- 
tematic comparative investigation of the merits of the dif- 
ferent solvers available for the quantum impurity model 
is needed. The relation between cluster geometry and 
molecule geometry should be more fully explored. Den- 
sity functional or GW methods may be better ways to 
treat the interactions not included in our DMFT calcu- 
lation. Also, alternative formulations of cluster methods 
better suited to finite systems may exist. 

Examination of more complicated systems, where the 
partitioning into correlated and "passive" subspaces and 
determination of which interactions to treat becomes 
more of an issue, is important. More generally, molecules 
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